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We report on the occurrence of small-scale vortices in simulations of the convective solar surface. Using an eigenanalysis of the 
velocity gradient tensor, we find the subset of high-vorticity regions in which the plasma is swirling. The swirling regions form 
an unsteady, tangled network of filaments in the turbulent downflow lanes. Near-surface vertical vortices are underdense and cause 
a local depression of the optical surface. They are potentially observable as bright points in the dark intergranular lanes. Vortex 
features typically exist for a few minutes, during which they are moved and twisted by the motion of the ambient plasma. The bigger 
vortices found in the simulations are possibly, but not necessarily, related to observations of granular- scale spiraling pathlines in "cork 
animations" or feature tracking. 
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1. Introduction 

Solar surface convection is characterized by upflows of hot 
plasma in Mm- scale granules and downflows of radiatively 
cooled plasma in a network of relatively narrow intergranular 
lanes (for a review, see |Nordlund et al.|2Q09| ). The turbulent na- 
ture of the downflows, together with angular momentum con- 
servation, suggests the formation of vortical flows at scales on 
the order of the width of the downflow lanes or smaller. In fact, 
vortex structures regularly appear in simulations of s olar sur- 
face convection with or without magnetic fiel ds (e.g., [Stein & 
Nordlund 1998; Vogler 2004; Muthsam et al.|2QT0| 
et al.,2011, Shelyag et al. 2011 
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There is observational evidence for photospheric whirl flows 
on larger scales between 5 Mm and 20 Mm (Bra ndt et al.|1988( 
[Simon & Weiss] 19971 [Attie et al.|2009 ). Granular- scale vortical 
motions associated with intergranular downflows ("sink holes") 
were found only more recently on the basis of high-resolution 
observations ( [Bonet et aL| 2008, 2010; Varg as Dominguez et al. 
[201 1| ). In such studies, the horizontal velocities are inferred ei- 
ther by local correlation tracking or by direct tracin g of indi- 
vidual features, such as magnetic bright points (e .g., Berger & 
[TltIel[T996l [Balmace da eTaL][20T0l ). [Steiner et al.] (|2010) found 
dark lanes moving into granules in high-resolution images ob- 
tained with the balloon-borne telescope SUNRISE. By compar- 
ison with numerical simulations, they identified these lanes as 
horizontally oriented vortex tubes, presumably originating from 
the overturning motion at the boundaries of the intergranular 
lanes (cf. Stein & Nqrdlundl998 ). Small-scale swirl events were 



also reported in chromospheric observations ( | Wedemey er-B ohm 
& Rouppe van der Voort|2009 ); their re lation to surface convec 



tion is unclear (cf. [Carlsson et al.|2010[ ). 

In this paper, we report on the detection of sub-intergranular 
scale vortices in simulations. We begin with a description of our 
simulations and the technique used for vortex identification in 
Sect.|2] In Sect. [3] we present a statistical analysis of vortices in 
our simulations, study the physical properties of selected speci- 
mens, and examine possible observational signatures. We sum- 
marize our findings in Sect. |4] 



2. Methods 

2.1. Simulations of solar convection 



We used the MURaM code ( |Vogler[[2003l [Vogler et al.][2Q05] ) to 

solve the MHD equations in the context of solar surface convec- 
tion, including radiative energy transfer and an equation of state 
which incorporates the eflTects of partial ionization. The simula- 
tions are performed on a Cartesian grid, using a 4th order cen- 
tral diflTerence scheme for the spatial discretization and a short- 
characteristics scheme for the radiative transfer. At the bottom 
of the domain, open boundary conditions allow for the in- and 



outflow of matter: d^v 



for outflows and v 



x,y 







for inflows, i.e., inflows are vertical. The pressure and the en- 
tropy of inflowing matter are chosen such that the total mass in 
the domain is approximately constant and the radiative energy 
loss is consistent with the solar luminosity. The top boundary 
is closed and the horizontal boundaries are periodic. The mag- 
netic field is assumed to be vertical at the top and bottom bound- 
aries: Bx,y = d^Bz = 0. For the radiative energy transfer, a gray 
approximation is assumed. The equation of state being used is 
appropriate for the conditions in the solar photosphere and the 
near-surface convection zone. 

The parameters of the simulations considered here are listed 
in Table [T] Runs C and H cover the same vertical range at in- 
termediate and high resolutions, respectively. Both simulations 
include weak, initially random magnetic fields and h ave previ- 
ously been used to study near- surface dynam o action ([Vogler & 
Schiissler"2007'; 'Pi etarila Graham et al.|20TQl [MoU et al.|2QTll 
Their high resolution allows us to study the fine details of small 



Table 1. Parameters of the simulations. The vertical range is 
given with respect to the height at which the horizontal mean 
temperature is ^7000 K. 



Run 


Box size [Mm^] 


Cell size [km"*] 


Vertical range [km] 


C 


4.86 X 4.86 X 1.4 


7.5 X 7.5 X 10 


-890 <z<510 


D 


12 X 12x6.144 


20.83 x20.83 x 16 


-5560 < z < 584 


H 


4.86 X 4.86 X 1.4 


4x4x4 


-884 <z<516 
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Figure 1. Snapshot of the swirling strength (green volume rendering) and the optical surface color-coded with vertical velocity 
(downflows in red and upflows in blue) in Run C. The size of the box shown is 4.8 x 4.8 x 1.4 Mm^. The optical surface is hidden 
in the lower right quadrant, uncovering the swirling structure in the subsurface layers. 



vortices. Run D does not include a magnetic field and uses a 
lower resolution, allowing us to probe the convection zone at 
greater depths. 

2.2. Vbrfex identification 

Despite there being a clear intuitive notion of what a vortex is 
(viz., the rotation of fluid parcels about a possibly moving axis), 
it is surprisingly difficult to come up with a formal definition. 
In fact, contemporary research in fluid mechanics is still con- 
cerned with finding unambiguous vortex identification schemes 
(e.g., |Haller|2005t [g)lar 2007 ). Given its frequent use in the lit- 
erature, it should be noted that a high vorticity oj = \V xv\is not 
a suflficient indicator for the presence of a vortex, because oj can 
also be high in shear flows without rotation. 

To detect vo rtices, we use the s o-called "swirling strength", 
or /lei, criterion ( Zhou et aLl|1999| ), which is based on an eige- 
nanalysis of the velocity gradient tensor Uij := (djVi). Vortices 
are defined to be regions where U has a pair of complex conju- 
gate eigenvalues. A large unsigned imaginary part implies a 
region of "strong swirling", i.e., (part of) a vortex. In the case 
of rigid rotation, the value of Ts := InjAd gives the revolution 
period of the rotating flow. Unlike many other methods for vor- 
tex detection, the criterion is applicable also in the case of 
compressible hydrodynamics ( Kolar 2009) . 

Depending on the signs of the real eigenvalues and the real 
part of the complex eigenvalues Acv, we discriminate four types 
of vortices (Haimes & Ken wright 1999) : 

1. /Ir < 0, /Icr < 0: spiraling inward, converging 

2. /Ir > 0, /Icr < 0: spiraling inward, diverging 

3. /Ir < 0, /Icr > 0: spiraling outward, converging 

4. /Ir > 0, /Icr > 0: spiraling outward, diverging 



Here, converging and diverging refer to the fluid motion along 
the direction of the eigenvector corresponding to A^, which is 
identified as the direction of the vortex. Note that this direction 
is not necessarily orthogonal to the "swirling plane" spanned by 
the real and imaginary parts of the complex eigenvectors. We 
distinguish vortices by their inclination angle l with respect to 
the vertical direction. 

The identification of vortices via A^ is Galilean invari- 
ant, because it only relies on derivatives of the velocity field. 
Streamlines, however, are not invariant under Galilean transfor- 
mations. Hence the streamlines in a particular frame of reference 
(e.g., the one which is at rest with respect to the computational 
box) do not always show circular orbits where vortices are de- 
tected. A strong vortex, however, is expected to produce vortical 
streamlines for a large range of observer velocities. 

Although the definition of the swirling strength is purely lo- 
cal (as is the definition of vorticity), it can be shown that for a 
volume of rigidly rotating fluid, the swirling strength is uniform 
over the entire interior of the volume. 



3. Results 

3.1. Occurrence of vortices 

In all our simulations, the regions of strong swirl form a net- 
work of highly tangled filaments. An illustrative example of this 
is shown in Fig.[T] Some of the filaments protrude from the opti- 
cal surface, forming variously shaped vortex tubes with diff'erent 
inclinations, sometimes bridging horizontally separated points in 
an arc. The optical surface is usually depressed at the footpoints 
of vertically emerging vortex tubes. 

The horizontal distribution of the swirls is not homogeneous, 
the granular upflows being mostly devoid of them. Near the op- 
tical surface, vortices with a large inclination with respect to the 
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veFticaf swirls (i < 20^} painted black horizontal gwirig (t > 70^) painted black 




Figure 2. Vertical velocities at the average height of the optical surface (T ^ 7000 K). Black spots indicate vortices with a swirling 
period Ts < 120 s and an inclination of less than 20° with respect to the vertical (horizontal) direction in the left-hand (right hand) 
panel. 




Figure 3. Horizontal area fraction of grid cells with swirling pe- 
riods below a given threshold Ts as a function of height in Run D. 
z = corresponds to the average height of the optical surface. 



vertical (l > 70°) are preferentially located at the edges of the in- 
tergranular downflow lanes, see Fig. [2 Those with a small incli- 
nation (l < 20°) are mostly inside the lanes where the downflow 
is strong. 

We find most of the strong vortical flows at a few hun- 
dred kilometers below the optical surface, see Fig. |3] With an 
increasing upper limit for the swirling period, the distribution 
declines less steeply with depth, i.e., most of the very deep 
swirls are slow. The structure at great depths is also filamentary. 
Horizontal swirls (dotted lines) are more numerous than verti- 
cal ones (dashed lines), consistent with isotropy. Near and above 
the optical surface, the distribution of inclination angles deviates 
from isotropy (see next Section). 



3.2. Statistical properties of vortex regions 

In Fig. |4j we present statistical properties of regions with large 
Aci. For these analyses, we combined data from five diff'erent 
snapshots of Run C with a separation of about 7 minutes each. 
All grid cells are treated equally in the statistics, regardless of 
whether or not they form part of a larger, contiguous swirling re- 
gion. In panels (b)-(h), we plot histograms normalized such as to 
give the fraction density with respect to the considered variable: 
the integral over a certain range of the variable gives the fraction 
of grid cells with values in that range. For every histogram we 
considered a 50 km (5 grid cells) vertical range about the indi- 
cated height, taking a total of -10^ grid cells into account. The 
individual panels of Fig. [4] are described in the following. 

(a) This plot shows, as a function of height, the fraction of 
space occupied by swirls with Ts < 120 s (compare Fig.jSj. The 
highest occurrence of these swirls is at z ~ -300 km (z = 
being the average height of the optical surface). Inwardly spiral- 
ing flows that diverge in the direction of the vortex (type 2, see 
Sect. 2.2 ) are prevailing everywhere; their mean fraction with re- 
spect to all types is 67% below z = and 75% above. The mean 
fraction of swirling flows of the outwardly spiraling, converging 
type (3) is 22% below z = -200 km and decreases to 7% above 
z = 100 km. Types (1) and (4) are insignificant. 

(b) This plot shows histograms of the swirling period Ts, nor- 
malized such that the integral over all Ts is one. The total frac- 
tion of grid cells with imaginary eigenvalues is about 53% at 
z = -90 km and z = 0, and 47% at z = 90 km. With decreasing 
height, the histogram is more peaked and the peak moves to- 
ward larger swirling strengths (smaller Ts). In all the other plots 
of Fig. [3[ only strong swirls with Ts < 120 s (i.e., to the left of 
the vertical line) are taken into account. 

(c) The histograms of the inclination angle l indicate that the 
swirls are preferentially (i.e., more than expected for an isotropic 
distribution) horizontal at z = and preferentially horizontal or 
vertical at z = 90 km, whereas at z = -90 km the distribution 
is consistent with an isotropic distribution (indicated by the gray 
dash-dotted line). 
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Figure 4. Vortex cell statistics for Run C. Except for panel (a), different heights are represented by brown {z - -90 km), green 
(z = 0) and blue lines (z = 90 km). Except for panel (b), only strong swirls with Ts < 120 s are taken into account. The dashed and 
dotted lines represent the subset of vortices with vertical {i < 20°) and horizontal {i > 70°) orientation, respectively. The fraction 
densities plotted in panels (b)-(f) are normalized histograms, the normalization being such that the integral over all values of the 
respective abscissa variable is one. Angle brackets (...) denote the horizontal mean of the respective variable in downflow regions 
(i.e., where < at the respective height). 
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Figure 5. 2D histograms of the swirling strength and various ob- 
servables (top panel: total magnetic field strength, middle panel: 
vertical magnetic field component, bottom panel: emerging in- 
tensity normalized by the overall horizontal mean) for Run C. 
All downflow regions in a ±25 km range about the average 
height of the optical surface are taken into account. Each his- 
togram is normalized such that an integral gives the fraction of 
volume with properties in the integrated range. 



(d) The gas pressure p in swirling regions is, on average, re- 
duced to 87% of the horizontal mean {p) in all downflow regions 
(v^ < 0) at the same height, consistent with a contribution by 
a dynamic pressure due to centrifugal forces. The gas pressure 
deficit does not depend much on height and inclination angle. 

(e) The swirling regions are significantly rarefied, the aver- 
age value of p/(p) being 84%. The rarefaction does not depend 
much on height and inclination angle. 

(f) The mean of the relative temperature T l{T) is 1.0 at all 
three considered heights. With increasing depth, the distribution 



becomes broader and the maximum moves towards lower tem- 
peratures. Most of the vertically oriented swirls (dashed lines) 
are slightly cooler than the horizontal ones (dotted lines), con- 
sistent with their predominant location inside cool intergranular 
lanes. 

(g) In horizontal swirls (dotted lines), the median and the 
mean of the vertical velocity are near zero, in agreement with 
them being found at the edges of the granules. In vertical swirls 
(dashed lines), the mean is roughly -3 km s"^ at all heights. 

(h) The mean of the horizontal velocity Vh is approximately 
3kms"\ irrespective of the height. Note that this velocity in- 
cludes the horizontal proper motion of the whole swirl. Vertical 
swirls tend to have slightly smaller horizontal velocities, pre- 
sumably because most of them (being at the center of intergran- 
ular lanes) are not subject to the bulk motion of the expanding 
granules. 

In Fig. |5j we present 2D histograms of the swirling strength 
and selected observables, taking into account all downflow re- 
gions near the optical surface. This restricts the analysis to 
the intergranular lanes. The top panel shows that strong mag- 
netic fields become increasingly rare with increasing swirling 
strength, there being no apparent correlation between the two. 
This is also true for both polarities of the vertical magnetic field 
alone. 

As is to be expected for the downflow lanes, the emerging in- 
tensity in vortex regions is preferentially smaller than the overall 
mean (bottom panel of Fig. [5]). There is no clear trend towards 
higher intensity with increasing swirling strength. This indicates 
that, although vortices are locally bright, the contrast is small 
compared to the brightness variations within intergranular lanes. 

3.3. Contiguous vortex features 

We define vortex features to be contiguous regions of grid cells 
with a large swirling strength. To achieve a clear separation, we 
consider only features above the (warped) optical surface and 
take a lower cutoff' for the swirling period than that used in the 
preceding section: z > z(t = 1) and Ts < 60s. In 6 indepen- 
dent snapshots of Run C we find 3806 features that satisfy these 
criteria (this number includes isolated grid cells). The number 
of features decreases rapidly with size, see panel (a) of Fig. [7 
It is roughly proportional to the inverse square of the occupiec 
volume (blue line). In the following, we consider only the 622 
features with at least 100 grid cells, corresponding to a minimum 
volume of 5.6 • 10"^ km^. 

Fig. [6] shows the features detected in one of the snapshots. 
On visual inspection, we find that the inclination angle l is in 
alignment with the longitudinal directions of the features, i.e., 
vertical (horizontal) features consist of swirls with small (large) 
inclination angles and more complicated shapes contain swirls 
with diff'erent inclinations. We use the median of the inclination 
angle l as a characteristic value for each feature. This median is 
denoted with t in the following. 

The histogram of t for all features is represented by the black 
line in panel (b) of Fig. [7] The red line represents the subset of 
features with a narrow distribution of angles, defined to be those 
for which 70% of all the constituent cells have an inclination l 
within ±10° about the feature's median t. The golden line repre- 
sents all other features. 

We divide the narrow subset (red) into two groups, a vertical 
one for which t is smaller than 45° (dashed lines) and a hori- 
zontal one for which t is larger than 45° (dotted lines). Density 
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Figure 6. Snapshot of Run C w ith t he warped optical surface in yellow and a volume rendering of the inclination angle l in se- 
lected vortex features (see Sect. 33 for a description of the selection criteria). Pink and green colors correspond to vertically and 
horizontally oriented vortices, respectively. The size of the box shown is 4.8 x 4.8 x 0.7 Mm^. 
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Figure 7. Statistical properties of contiguous vortex features above the optical surface in Run C. The horizontal number density 
plotted in panel (a) is the number of features per horizontal area. In panels (b)-(f), features are classified by the median inclination 
t and the width of the distribution about the median. Vertical (horizontal) features are deemed to be those with narrow distributions 
and small (large) r, they are represented by dashed (dotted) lines in panels (c)-(f). The plotted fraction densities are histograms 
normalized such that the integral over all values of the respective abscissa variable is one. In the abscissa labels, angle brackets (. . .) 
denote the (height-dependent) horizontal mean of the respective quantity in downflow regions and a bar denotes the mean within a 
particular feature. 
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Figure 8. A strong vertical vortex penetrating the optical surface in Run H. The nine panels at the top show, at various stages of 
development, the velocity (left panels) and gas pressure {middle panels) at the average height of the optical surface (z = 0) together 
with a bolometric intensity map {right panels). The line plots at the bottom left represent one-dimensional cuts through the vortex at 
the last instant {t = 126 s). The 3D visualization at the bottom right shows, from two different viewing angles, selected streamlines 
color coded by velocity and the warped optical surface in yellow; the black rectangle indicates the z = plane. 
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Figure 9. Formation of a strong vertical vortex in Run C. The plots display the swirling strength (green volume rendering) and the 
optical surface (yellow) at three different times (labels are in minutes). The size of the box shown is 1.1 x 1.1 xO.8 Mm^. 




Figure 10. Rise and fall of a vortex arc in Run C. The plots display the swirling strength (green volume rendering) and the optical 
surface (yellow) at three different times (labels are in minutes). The size of the box shown is 1.5 x 1.5 x 0.8 Mm^. 



functions of the gas pressure, density, temperature and vertical 
velocity are plotted in panels (c)-(f); each quantity is divided by 
the horizontal mean in downflow regions at the respective height 
and then the mean within the respective feature is considered. 
The results are consistent with the "point-by-point" statistics 
presented in Sect. 32^ (compare with Fig. |4]): gas pressure and 
density are significantly reduced with respect to the horizontal 
mean, the differences between the horizontal and vertical sub- 
set being qualitatively similar. The vertical velocity is strongly 
negative (downward) in almost all of the vertical features. 

Lifetimes Through visual inspection, we estimate the lifetimes 
of 17 features in Run C. We define the lifetime to be the time 
interval in which a feature is visible as a distinct entity above the 
optical surface. The resulting mean lifetime is 3:30 minutes, the 
standard deviation being 1:40 minutes. 

3.4. Individual examples and evolution 
3.4.1. Vertical vortices 

Fig. [8] displays an example of a vertically oriented vortex which 
penetrates the optical surface. Starting out as a somewhat dis- 
torted vortical flow (in the box's frame of reference), it evolves 
into a roundish vortex with a diameter of approximately 80 km. 
The rotation in the horizontal plane is nearly rigid at its core 
(innermost 60 km) and the vertical velocity reaches a local max- 
imum at the center of the vortex. Both gas pressure and density 
are lowered to as little as 36% of the respective horizontal mean 
at z = and are significantly distinct from the local background. 
The optical surface is depressed inside the vortex, the bottom of 
the depression being 110 km below z = 0. At the site of the vor- 



tex, the bolometric intensity is locally increased by up to 24% 
with respect to the mean. The vortex thus appears as a bright 
spot within the dark intergranular lane. 

Fig. [9] illustrates the formation of a vertical vortex. The fea- 
ture is indicated by a black arrow in the three snapshots shown. 
It starts as a conglomeration of weakly swirling structures that 
apparently merge into a single distinct feature. The swirling 
strength increases with time. At all stages, the optical surface 
is depressed in the vortex, the largest depression being on the 
order of -100 km. Above the optical surface, the feature has a 
length of -350 km. 

3.4.2. Vortex arcs 



Fig. 10 displays the formation of a vortex arc. When the feature 
emerges through the optical surface, it is rather weak and hard to 
distinguish from the background. As it rises to its peak height, it 
gains in swirling strength. The horizontal distance between the 
two footpoints on the optical surface is -780 km at this stage, 
the height above the surface is -230 km. At last, the feature de- 
velops a kink near its crest and collapses. In total, it lives for 
approximately 4 minutes. 

Fig. [TT] depicts a section from a snapshot of Run H which 
contains a big, slowly evolving vortex arc that hovers high above 
the optical surface. Next to it, on the right-hand side in the plots, 
is a much smaller and rapidly moving vortex arc. The two arcs 
are indicated with red arrows in panel (a). In the case of the big 
arc, the separation of the footpoints of the streamlines on the op- 
tical surface is -1.4 Mm, the peak height above the optical sur- 
face is -330 km. The vortical fluid motion, indicated by the large 
swirling strength in panel (a), is well visible in the streamlines 
plotted in panels (b) and (c). The purple and turquoise colors in- 
dicate the horizontal direction of the fluid motion: material flows 
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from outside towards the big arc's crest and down on both of its 
legs. At the part of the big arc where the streamlines are purple 
on one side and turquoise on the other, it is subject to longitudi- 
nal shear. 

The small vortex arc has a horizontal extent of -310 km and 
a peak height of -210 km above the optical surface. Moving 
rapidly with respect to the computational frame of reference, 
it does not show vortical streamlines in panels (b) and (c) of 
Fig.[TT] Its evolution is depicted in Fig. [12] In 75 s, it propagates 
-200 km horizontally and rises -100 km vertically. The corre- 
sponding velocities are 2.7 km s"^ and 1.3kms"\ respectively. 
The streamlines show that it corresponds to the front of an up- 
ward flow. In the upper, horizontal part of the arc, the mean den- 
sity is reduced to 87% of the overall horizontal mean. A crude es- 
timate of the buo yant rise velocity limited by aerodynamic drag 
(Eq. 6 in |Parker|l975) yields -1.6 km s ^ While this is consis- 
tent with the measured vertical velocity, it is not clear whether 
buo yancy plays a n important active role in the vertical motion 
(cf. |Arendt|1993] ). 

3.5. Relationship with larger-scale vortical motion 

As described in the introduction, observational evidence for vor- 
tical motion on the solar surface typically refers to larger scales 
than the strong, small-scale vortices studied here, although the 
general characteristics are similar (e.g., the association with 
downflows). The reported mean lifetimes of granular-scale vor- 
tices of 5-8 minutes ( Bonet et al.|2008[ |2010) are not drastically 
difl'erent from the value of 3.5 min estimated above. However, 
because their rotation periods are difl'erent, our vortices make 
about 2 revolutions during their lifetimes, while the observed 
vortices can be followed for only a fraction of one rotation (for 
instance, about 25% of a period in the case of Bo net et al.|2010 ). 
It is conceivable that the observed vortical motions represent the 
peripheral parts of the much stronger small-scale vortex cores 
that show up in the simulations but are too small to be observed 
directly. The outer vortex parts would be much more strongly af- 
fected by the evolving granulation pattern and thus be detectable 
only for a fraction of a rotation period. 

To see whether the vortices studied in this paper would be 
detectable in observational data through feature tracking tech- 
niques, we consider horizontal pseudo pathlines (trajectories of 
fluid elements) in Fig. 13 The pathlines are determined from 10 
snapshots of the horizontal velocity field at the average height of 
the optical surface with a temporal spacing of ^^18 s. They are 
"pseudo" because the vertical velocity is ignored. 

Some pathlines are spiraling in towards a region of high 
swirling strength, see the feature at (x,y) ^ (3.7, 1.9) Mm in 
the left-hand panel of Fig. 13_ In this case, there is a clear as- 
sociation with the small-scale vortices presented in this paper. 
However, we also find pathlines which are curved but for which 
the association with an actual vortex is not clear, for instance in 



the blown-up region shown in Fig. 13 



4. Summary and discussion 

We have investigated vortical fluid motions in simulations of 
near-surface solar convection by calculating the eigenvalues 
and eigenvectors of the velocity gradient tensor field. Complex 
eigenvalues with a large imaginary part indicate regions of 
strong swirling. They are found predominantly in and near the 
intergranular lanes, where cooled fluid is sinking down in a tur- 
bulent fashion. The swirling regions form an unsteady network 



of highly tangled filaments, some of which protrude above the 
optical surface. 

Near the optical surface, vertically oriented swirls are pref- 
erentially located in the interior of intergranular lanes, where the 
downflow is strong. Horizontal swirls, on the other hand, are 
predominantly located at the edges of the granules, where ver- 
tical motion is mostly absent. The 3D structure of contiguous 
features above the optical surface is manifold, but often in the 
form of bent and arc-shaped filaments. These type of structures 
have previously been seen in independent numerical simulations 
with diff'erent codes, notably Muthsam et al. ( 2010) and Stei n &| 
Nordlund (1998). The swirling direction (rotation axis) is typi- 
cally aligned with the longitudinal direction(s) in a contiguous 
feature. 

Rotary motion causes dynamic pressure by centrifugal 
forces. One may, therefore, expect that vortex tubes are under- 
dense as a result of pressure equilibrium and thermalization with 
the environment, analogous to magnetic flux tubes. We find vor- 
tex features above the optical surface to be underdense by -15% 
with respect to the horizontal mean in downflow regions. A sig- 
nificant depression of the optical surface is common where it 
intersects with vortices. Having diameters of -80 km, the largest 
vertical vortex features would be visible as bright spots with a 
size of -0.1 arcseconds in observations of the Sun. 

Most vortex features are highly unsteady, being moved 
and/or twisted by turbulent fluid motions. Some of them rise up- 
war ds. Although, in principle, buoyancy coul d contribute to this 
rise ^Parker! 1991 1 but see also |Arendt||1993| ) we find no unam- 
biguous evidence for the operation of buoyancy forces. After a 
typical lifetime of a few minutes, the features usually descend 
below the optical surface, where they dissolve. 

In general, the vortices presented here are smaller than those 
found in observations so far. They may represent peripheral parts 
of strong, small-scale vortices. However, arc-like motions of sur- 
face features do not necessarily imply the presence of a strong 
central vortex. The pattern of the fluid flow is changing consid- 
erably at time scales of only a few minutes, which is reflected by 
the lifetimes of the small-scale vortices. We therefore reckon that 
granular- scale or bigger vortical flows are less numerous than 
small-scale vortices. 
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Figure 12. Evolution of the small vortex arc which is shown in a 
blow-up in Fig.fTT] The isosurfaces show the arc at four different 
times with a separation of 25 seconds. For each stage, selected 
velocity streamlines are plotted in the corresponding color. 



Figure 11. Snapshot of Run H with vortex arcs. All plots dis- 
play the warped r = I surface in gray color. Panel (a) shows 
the swirling strength (green volume rendering), red arrows indi- 
cate the features discussed in the text. Panels (b) and (c) show 
selected streamlines and the vertical velocity (red/blue volume 
rendering). The streamlines are colored according to the hori- 
zontal velocity in the approximate direction of the big arc. The 
size of the box shown is 1.3 x 1.4 x 0.7 Mm^. The blow-up in 
panel (a) contains the feature presented in Fig. [12 
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